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ANNUAL  REPORT 


FEBRUARY  1993  -  JANUARY  1994 
Early  Detection  of  Breast  Cancer  on  Mammograms:  Perceptual  Feedback, 
Computer  Processed  Images  and  Ultrasound 

Grant  #  DAMD17-93-J-3014 
Principal  Investigator:  Peter  Bloch,Ph.D.: 

The  treatment  of  small  breast  lesions  (less  than  1cm)  with  lumpectomy  and 
definitive  radiation  therapy  results  in  excellent  clinical  results  with  10  year  disease  free 
survival  of  better  than  95%.  The  project  is  to  improve  the  early  detection  of  breast  cancer 
when  the  lesions  are  too  small  to  be  palpable  when  the  treatment  is  most  likely  to  be 
effective  in  sterilizing  the  disease.  Three  approaches  are  being  investigated  as  a  method  for 
improving  the  detection  of  small  breast  malignancies.  These  include;  using  (1)  Perceptual 
Feedback  as  an  Aid  to  the  Early  Detection  of  Cancer  on  Mammograms,  (2)  Digital  Image 
Processing  of  Mammograms,  and  (3)  High  Resolution  Ultrasound  Mammography.  The 
progress  achieved  during  the  year  using  these  three  approaches  are  summarized  in  this 
annual  report. 

Project  lA:  Perceptual  Feedback  as  an  Aid  to  the  Early  Detection  of  Cancer 
on  Mammograms. 

Co-Investigator:  Harold  L.  Kundel,  M.D. 

INTRODUCTION:  RESEARCH  OBJECTIVES 

About  15%  of  tumors  that  are  actually  visible  in  images  are  missed  on  the  initial 
reading  of  screening  mammograms.  These  misses  are  considered  to  be  perceptual  errors. 
Eye-position  recordings  during  the  search  for  lung  tumors  on  chest  radiographs  have 
shown  that  about  15%  of  perceptual  errors  are  due  to  incomplete  scanning  of  the  image 
with  the  eyes,  15%  are  due  to  a  recognition  mechanism,  and  70%  are  due  to  incorrect 
decisions.  Eye-position  feedback  which  consists  of  showing  where  the  gaze  has  been 
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concentrated  has  also  been  shown  to  decrease  the  number  of  incorrect  decisions.  This 
project  will  use  the  same  techniques  that  have  been  used  for  chest  radiographs  on 
mammograms.  The  following  methodology  is  being  followed.  (1)  Develop  a  computer 
display  for  mammograms  that  matches  the  grayscale  properties  of  the  display  to  the  contrast 
discrimination  function  of  the  observers  eye  and  that  displays  multiple  images  for  within 
breast  and  between  breast  comparison.  (2)  Interface  a  head-eye-position  recording  system 
to  the  display.  (3)  Develop  and  digitize  a  set  of  mammograms  with  barely  visible  visual 
cues  for  tumor.  (4)  In  a  case-control  study,  compare  detection  performance  with  and 
without  eye-position  feedback. 

BODY:  PROGRESS  REPORT 

The  computer  display  has  been  developed  using  a  SUN  SPARC  10  platform.  Dome 
video  controller  boards  and  Tektronix  monitors.  The  monitors  have  a  maximum  luminance 
of  80  fl,  a  2200  X  1800  X  8  matrix  and  a  191  micron  (7.5  mil)  scanning  spot.  The  major 
advantage  of  the  monitors  is  the  luminance  and  the  scanning  spot  uniformity.  A  display 
format  that  shows  two  related  views  and  a  set  of  three  icon  pairs  on  one  monitor  has  been 
selected.  Two  monitors  can  then  be  used  to  show  either  one  two-view  study  (R  and  L 
sides)  or  to  compare  recent  and  old  images  for  one  side.  The  icons  are  small  versions  of 
other  views  and  other  studies.  The  icons  are  used  by  the  operator  to  select  other  studies  for 
comparison  and  review. 


Monitor  1  Monitor  2 
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The  head  position  sensor  obtained  from  the  Applied  Science  Laboratories  (ASL) 
has  been  interfaced  successfully  to  an  existing  eye-position  monitor.  Eye-scanpaths  have 
been  recorded  to  test  the  system  and  to  provide  preliminary  data  for  the  analysis  system.  An 
example  of  a  20  second  segment  of  an  eye-position  pattern  is  shown  in  the  figure  below.  A 
diagram  of  the  mammogram  is  on  the  left  side.  It  illustrates  a  mass  in  the  left  breast.  The 
eye  scanpath  of  an  observer  is  on  the  right  There  is  considerably  more  comparison  activity 
than  has  been  observed  for  chest  images  and  as  with  chest  images  there  is  focal  attention  on 
suspected  abnormalities  as  shown  by  the  concentration  of  fixations  on  the  left  breast  mass. 


Mammographic  Display  with 
a  Mass  in  the  Left  Breast 


CONCLUSIONS: 

The  research  program  is  progressing  as  planned.  During  the  next  grant  year,  we 
will  concentrate  on  (1)  displaying  mammograms,  (2)  obtaining  a  test  set  of  digitized 
mammograms  and  (3)  measuring  eye-position  parameters  such  as  fixation-cluster 
distribution  and  cluster  dwell  time.  Displaying  mammograms  effectively  is  a  non-trivial 
matter  and  requires  great  attention  to  detail.  The  input-output  transfer  characteristic  (lOTC) 
of  the  monitors  needs  to  be  calibrated  carefully  and  matched  to  the  observers  contrast 
sensitivity  function  at  the  ambient  luminance  level.  It  may  be  necessary  to  display  the  same 
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breast  images  using  different  lOTCs  in  order  to  take  full  advantage  of  the  dynamic  range 
(12-14  bits)  of  the  digital  image.  Unless  the  mammograms  are  properly  displayed,  the  use 
of  perceptual  feedback  will  be  superfluous.  There  is  still  some  development  worit  that 
needs  to  be  done  on  this  problem.  Dr.  Orel  is  already  working  on  the  collection  of  images 
but  they  still  require  digitization.  There  is  very  little  development  necessary  on  the  eye- 
position  analysis  programs.  The  major  task  will  be  to  determine  if  as  in  the  case  of  lung 
tumors  and  fractures  missed  lesions  receive  extra  focal  visual  attention.  This  information 
will  then  be  used  to  set  the  thresholds  in  the  perceptual  feedback  programs. 

Project  IB:  Digital  Image  Processing  of  Mammography  Films 
Co-Investigator:  E.  Loren  Buhle,  Jr.,  Ph.D. 

INTRODUCTION:  RESEARCH  OBJECTIVES 

The  treatment  of  breast  cancer  has  a  much  higher  success  rate  if  the  carcinoma  is 
detected  at  early  stage.  Mammography  is  the  only  diagnostic  procedure  with  a  proven 
capability  for  detecting  early-stage,  clinically  occult  breast  cancers  [  ACS,  1983]. 
Diagnostic  information  can  come  from  both  stromal  or  parenchymal  patterns,  contained  in 
the  low  spatial  resolution  domain  and  clusters  of  microcalcifications,  found  in  the  high 
spatial  resolution  domain. 

Many  investigators  have  attempted  to  correlate  parenchymal  patterns  with  breast 
cancer  [Wolfe,  1976,  1978;  Wellings  et  al.,  1978;  Saftlas  et  al.,  1987].  The  major 
components  of  the  parenchymal  patterns  are  the  amounts  of  nodular  and  homogeneous 
densities  seen  on  the  mammogram  [Brisson  et  al.,  1982,  1984].  While  controversial, 
breast  cancer  risk  has  been  suggested  to  increase  progressively  with  the  extent  of  these 
densities  on  the  mammogram  {Brisson  et  al.,  1982, 1984;  Wolfe  et  al.,  1987],  although 
this  association  seems  stronger  for  nodular  densities  compared  to  homogeneous  density 
[Brisson  et  al.,  1982, 1984].  The  parenchymal  pattern  and  percentages  of  the  breast 
showing  densities  are  also  related  to  other  factors,  such  as  obesity  and  patient  age. 
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Approximately  30-50%  of  the  breast  carcinomas  reveal  microcalcifications  on 
mammograms,  and  between  60-80%  of  the  breast  carcinomas  reveal  microcalcificaticms 
upon  microscopic  examination  [Murphy  et  al.,  1978].  Therefore,  any  increase  in  the 
detection  of  microcalcifications  has  a  strong  potential  to  improving  the  efficacy  of 
mammography  for  detecting  ea.rly  breast  cancer.  The  efficiency  and  effectiveness  of  a 
mammography  screening  program  may  be  greatly  increased  with  computer  image 
processing  can  be  employed  for  the  detection  of  malignancy. 

Several  investigators  have  used  various  computational  schemes  to  analyze  breast 
abnormalities  from  the  low  resolution  regional  density  pattern  component  of  mammograms. 
Winsberg  et  al.  [1967]  compared  density  patterns  within  individual  breasts  and  between  the 
left  and  right  breasts.  Ackerman  et  al.  [1972]  used  Bayesian  probabilities  techniques  on  the 
regional  tissue  texture  found  in  xeromammograms  to  classify  benign  and  malignant  lesions. 
Ackerman  optimized  the  feature  space  in  order  to  minimize  the  misclassification  rate  for  a 
two  class  problem  by  trial  and  error.  This  work  spurred  further  investigations  by  Kimme 
et  al.  [1979],  Semmlow  et  al.  [1980],  and  Magnin  et  al.  [1986]  to  perform  statistical 
analysis  of  textural  features  of  breast  images  and  comparing  corresponding  features  in  the 
left  and  right  breasts.  This  work  used  small  data  sets  and  resulted  in  numerous  false 
positives. 

The  effectiveness  of  film-screen  mammography  for  screening  malignant  disease 
may  be  significantly  increased  by  computer-aided  detection  of  microcalcifications.  Recent 
studies  at  the  University  of  Chicago  [Zhang  et  al.,  1983;  Nishikawa  et  al.,  1993;  Ema  et 
al.,  1993]  reported  85%  true  microcalcification  cluster  detection  with  computer-aided 
detection  with  a  false-positive  detection  ranging  from  0.4-0.7  per  image,  using 
mammograms  digitized  at  l(X)-175n  at  10-bits  of  dynamic  range.  Digitizing  the 
mammogram  at  42[i  and  16-bit  dynamic  range  would  preserve  the  spatial  resolution  and 
latitude  already  present  in  the  film-screen  combination.  This  potential  for  increased  true¬ 
positive  and  decreased  false-positive  detection  still  needs  to  be  determined.  The  importance 
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of  maintaining  the  wide  latitude  of  the  original  film-screen  mammogram  in  the  digitized 
image  will  be  tested  by  evaluating  the  computer  aided  detection  of  microcalcifications  as  a 
function  of  the  regional  optical  density  of  the  mammogram. 

The  target-film-screen  combination  presently  employed  in  mammography  results  in 
a  spatial  resolution  of  approximately  20  line-pairs/mm  or  50  microns  [Lawrence,  1977]. 

To  preserve  this  resolution  in  a  digitized  image,  the  sampling  optimally  should  be 
approximately  25^^  Murphy  et  al.  [1978]  detected  many  more  microcalcifications  from 
histopathological  sections  than  what  were  seen  in  the  ccnresponding  mammograms.  Many 
of  the  studies  digitized  with  a  dynamic  intensity  range  of  255  gray  levels.  The  ability  to 
digitize  the  full  dynamic  range  of  the  mammogram  and  allow  a  user-specified  gray  level 
window  (e.g.  as  used  in  visualizing  CT  sections)  would  facilitate  the  diagnostic  evaluation 
of  mammograms  containing  regions  of  high  and  low  densities. 

The  hypothesis  to  be  tested  in  this  project  are;  Computer  processing  of 
mammograms  can  increase  the  observer  identification  of  nucrocalcifications  in  breast 
tissues.  Characterization  of  the  breast  parenchyma,  correlated  with  the  location  of 
microcalcification  clusters  can  be  used  to  identify  the  location  of  suspected  breast 
malignancies.  Combining  the  low  resolution  information  regarding  breast  parenchyma 
with  the  high  resolution  features  from  microcalcification  can  be  useful  in  detecting  and 
distinguishing  benign  and  malignant  lesions  of  the  breast 
BODY:  PROGRESS  DURING  HRST  YEAR 

The  first  year  was  devoted  to  examining  commercially  available  digitizers  capable  of 
high  spatial  resolution  and  a  wide  dynamic  range.  A  test  suite  of  phantom  test  objects  was 
imaged  using  conventional  mammogram  techniques  and  the  modulation  transfer  function 
(MTF)  of  the  x-ray  target-screen-film  combination  and  the  film  digitizer  was  evaluated. 

^The  Nyquist  sampling  criterion  asserts  that  a  continuous  waveform,  when  sampled  at  a  frequency  greater 
than  twice  the  maximum  frequency  of  the  object  being  measured,  can  be  reconstructed  completely  from  the 
sampled  waveform.  Conversely,  if  a  continuous  waveform  is  sampled  at  a  frequency  lower  than  twice  its 
maximum  frequency  component,  a  distrotion  termed  aliasing  occurs.  This  is  often  called  undersampling. 
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The  MTF  and  the  noise  properties  of  the  final  imaging  system  will  be  employed  to  design 
an  algorithm  to  restore  some  of  the  high  spatial  frequency  components  lost  during  image 
acquisition.  The  MTF  of  the  mammogram-digitizer  combination  was  measured  at  the  100, 
50,  and  20\i,  using  a  high  resolution  densitometer  present  in  the  Department  of  Cell  and 
Development  Biology  at  the  University  of  Pennsylvania. 

The  transfer  of  imaging  information  occurs  in  two  stages:  exposure  of  the  film- 
screen  combination  of  the  mammographic  acquisition  system  and  the  digitization  of  these 
films.  This  combined  transfer  function  of  the  digitized  film  was  measured  for  image  pixel 
sizes  ranging  from  100-10  microns.^  The  MTF  was  determined  for  digital  image 
acquistion  employing  different  pixel  sizes  by  analyzing  a  step  function  present  in  the  ACR 
mammography  test  phantom.  The  derivative  of  the  step  function,which  is  the  line-spread- 
function  (LSF)  was  accurately  represented  by  a  Gaussian  function.  The  modulus  of  the 
Fourier  Transform  of  the  LSF  is  the  MTF  of  the  system.  The  results  of  these  studies 
indicated  that  the  resolution  limit  (assumed  to  be  at  the  10%  level  of  the  MTF)  was  2.98 
line-pairs/mm  (335|i)  for  the  film  scanned  for  lOOfi  pixels,  10.62  line-pairs/mm  (94^)  for 
50)1  pixels  and  1 1.52  line-pairs/mm  (86)1)  for  20)i  pixels.  While  this  densitometer  would 
not  be  used  for  routine  digitization  of  mammograms,  it  was  clear  the  prospective  digitizer 
should  be  capable  of  sampling  at  a  50)1  rate  or  bener.  It  was  noted  that  dynamic  range  of 
the  mammograms  routinely  exceeded  8-bits,  ranging  from  !0  to  12  bits  contrast. 

A  commercial  digitizer  recently  became  avmlable  that  is  suitable  to  this  criteria.  We 
are  currently  evaluating  the  DBA  Systems,  Inc.  scanner  and  will  revise  the  resolution 
measurements  obtained  above.  The  noise  characteristics  will  also  be  evaluated  at  this  time. 

Preliminary  information  from  the  DBA  Systems,  Inc.  scanner^  was  also  used  to 
being  the  acquisition  of  computer  hardware  and  the  development  of  software  tools  in 
archiving  and  processing  this  massive  amount  of  information.  A  Digital  Equipment 

^Using  a  Perkin-EImer  flatbed  scanner 
^DBA  Systems,  Inc.,  Melbourne,  Florida 
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Corporation  3000-500  workstation  was  purchased,  with  the  appropriate  network  interfaces 
to  allow  other  researchers  direct  access  to  the  scanned  images.  Because  of  the  rapid 
development  of  the  hierarchical  storage  memory,  purchase  of  the  extensive  storage 
hardware  is  being  deferred  until  early  1994.  The  software  composing  the  patient's  archive 
is  based  on  the  treatment  planning  archive'^  discussed  elsewhere  [Buhle,  1993]  and  was 
extended  to  take  advantage  of  the  characteristics  of  the  DBA  image  information.  Using  the 
IDL  programming  language^,  software  tools  were  developed  to  selectively  enhance  the 
high  spatial  frequencies  information  in  the  digitized  mammogram,  while  maintaining  an 
optimal  signal-to-noise  ratio  and  thus  assist  in  the  computer  aided  detection  of 
microcalcifications.  Morphological  operators  and  region-growing  software  will  be 
developed  to  quantify  the  irregularly  shaped  microcalcification  clusters.  Software  tools  are 
also  under  development  to  describe  and  quantify  the  density,  texture  and  overall  disposition 
of  the  breast  parenchyma. 

CONCLUSION: 

The  tools  developed  in  the  first  year  will  be  applied  to  a  mixture  of  approximately 
80-100  mammogram  studies  containing  both  benign  and  malignant  pathology.  The  low 
resolution  information  will  be  examined  in  an  attempt  to  describe  and  quantify  a  normal 
parenchymal  pattern  of  the  breast,  beyond  the  classical  categories  suggested  by  Wolfe 
[1976,  1978].  There  is  no  current  description  of  the  normal  parenchymal  pattern. 

The  high  resolution  information  contained  in  the  digitized  mammogram,  band- 
limited  to  maximize  the  signal-to-noise  ratio,  will  be  examined  for  clusters  of 
microcalcifications.  Region  growing  software  routine  [Gupta  et  al.,  1987]  will  be 


'^The  patient  archive  is  based  on  the  public  domanin  software  [Rew  et  al.,  1990]  and  allows  the  construction 
of  centralized,  machine  independent  image  and  non-image  information.  The  self-describing  information 
present  in  this  archive  is  readily  assessible  via  C,  Fortran,  IDL,  and  C++  programs  running  on  a 
hetereogeneous  computing  enviroi.ment,  without  specific  knowledge  of  the  data  format  or  content  of  the 
archive.  The  self-describing  nature  and  extensible  nature  of  this  archive  permits  degrees  of  flexibility 
reflecting  the  unique  nature  of  each  patient's  study.  A  number  of  two-dimensional  mammograms,  three 
dimensional  image  volumes  (e.g.  CT  or  MRI),  geometric  contours,  and  non-image  infomration  can  be 
stored  in  each  archive  and  be  assessed  by  other  programs  running  in  different  computing  environments. 

^IDL  -  Interactive  Data  Language,  from  Research  Systems,  Inc.,  Boulder,  Colorado 
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implemented  to  describe  the  regularity  of  the  outer  border  of  these  clusters.  The  overall 
locations  of  these  microcalcifications  clusters,  relative  to  the  breast  parenchyma  and  other 
descriptive  elements  of  the  breast  (e.g.  nipple,  chest  wall),  venous  diameter,  displacement 
of  parenchyma  by  transparent  lesions  will  also  be  examined. 
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Project  1C:  High  Resolution  Ultrasound  Mammography 
Co-Investiagtor:  Bernard  Steinberg,  Ph.D. 

INTRODUCTION 

The  objective  of  the  ultrasound  echoscanning  mammography  research  is  detection 
and  classification  of  breast  lesions  the  order  of  2  mm  in  size  and  differentiation  of  such 
small  tumors  from  cysts.  The  primary  obstacle  is  poor  image  quality  (artifacts,  false 
targets)  caused  by  distortion  to  the  wavefront  when  the  transducer  is  made  large  enough  (5- 
10  cm)  to  achieve  the  desired  lateral  resolution  (0.2  mm)  and  sidelobe  level  (~  -60  dB)  at 
chest  wall.  In  a  prior  experimental  study  we  obtained  a  massive  database  on  wavefront 
distortion  of  acoustic  waves  as  they  passed  through  in  vivo  breast.  Analysis  of  these  data 
was  the  first  essential  task  under  this  program,  to  be  followed  by  algorithmic  development 
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of  signal  processing  methods  to  counter  wavefront  distortion,  and  the  development  of 
sensitive  experimental  facilities  to  measure  the  effectiveness  of  the  procedures. 

BODY;  SUMMARY  OF  TASKS 

1 .  In  vivo  data  analysis.  450  MB  of  data  had  been  measured  earlier  on  100  breasts, 
largely  of  mammography  patients  at  the  Hospital  of  the  University  of  Pennsylvania  (HUP). 
The  breasts  were  categorized  according  to  premenopausal  dense,  premenopausal  fatty,  and 
postmenopausal.  An  extensive  analysis  was  performed  this  year  resulting  in  several 
journal  articles,  of  which  three  have  been  published  and  four  accepted  subject  to  revision, 
and  five  papers  given  at  professional  society  meetings.  Details  of  the  analysis  are  found  in 
the  Appendix. 

2 .  Evaluation  of  existing  wavefront  distortion  compensation  algorithms.  A  large 
number  of  algorithms  have  been  published  by  our  laboratory  and  others.  All  or  nearly  all 
can  be  shown  to  be  phase  correction  algorithms,  that  is,  they  compensate  for  the  distorting 
to  the  phase  of  the  received  waveform  but  not  the  amplitude  or  modulus.  We  evaluated 
phase  aberration  correction  on  in  vivo  data.  The  results  are  helpful  but  inadequate;  these 
algorithms,  which  all  operate  upon  the  transverse  spatial  correlation  function  of  the 
distorted  wavefront,  reconstruct  a  main  beam  reasonably  well,  but  do  not  result  in  the  low 
sidelobes  required  for  successful  differentiation  of  tumors  from  cysts.  Thus,  this  genre  of 
algorithms,  which  is  the  largest  class  of  existing  algorithms,  are  inadequate. 

3 .  Development  of  newer  andlor  stronger  algorithms.  The  data  analysis  produced 
evidence  that  a  major  distortion  inducing  phenomenon  within  the  breast  is  refraction,  in 
contrast  to  weak  scattering  which  hitherto  had  been  thought  to  be  the  dominant  cause  of  the 
trouble.  Weak  scattering  responds  very  favorably  to  phase  aberration  correction. 
Refraction  does  not.  We  believe  that  much  of  the  refraction-induced  wavefront  distortion 
occurs  at  the  fat-glandular  tissue  interface.  We  are  working  on  three  algorithmic 
procedures  that  address  this  problem. 
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4.  Acquisition  of2-D  data.  The  wavefront  distortion  sources,  whether  scattering  or 
refraction,  disperse  energy  volumetrically  and  not  only  in  the  plane  of  the  transducer.  Our 
data  analysis  showed  that  the  algorithms  required  would  necessarily  be  2-D  algorithms. 
The  HUP  measurements  upon  which  these  observations  were  made  were  in  one  dimension 
only.  To  insure  that  the  conclusions  drawn  from  the  1-D  data  are  sound,  it  is  necessary  to 
acquire  2-D  wavefront  data.  It  is  also  necessary  to  have  such  data  to  test  the  algorithms  that 
will  develop  from  the  data  analyses.  For  this  reason  we  developed  a  collaboration  with 
Professor  Robert  Waag  at  the  University  of  Rochester  who  was  set  up  to  make  2-D 
wavefront  measurements  of  acoustic  waves  passing  through  in  vitro  breast  specimen. 
Work  of  this  nature  was  originally  scheduled  for  about  a  year  from  now,  after  we  had 
developed  our  own  experimental  facility.  For  the  reason  given  above,  however,  we 
decided  to  accelerate  the  process.  We  took  data  in  Professor  Waag's  laboratory  this  past 
summer  and  the  data  analysis  is  under  way. 

5  .  2-D  transducer  array  analysis.  Because  of  the  conclusion  reported  above  that  the 

algorithms  required  would  necessarily  be  2-D  algorithms,  we  have  begun  an  analytical  and 
computer  simulation  study  of  large,  2-D  transducer  arrays  (up  to  10  cm),  highly  thinned(or 
sparse)  for  practicality  and  cost  control.  Preliminary  results  favor  the  highly  thinned, 
deterministic,  multiple  ring  array  over  the  equally  highly  thinned  random  array. 

6.  Experimental  facility.  Development  of  a  3-D  synthetic  aperture  water  tank 
experimental  facility  is  under  way.  It  will  serve  two  functions:  Precision  measurement  of 
radiation  patterns  of  large  (up  to  10  cm),  2-D,  highly  thinned  transducer  arrays;  and  echo 
tests  of  deaberration  algorithms  on  phantoms.  We  expect  to  have  the  assembly  completed 
within  6-9  months.  We  will  then  evaluate  the  system  to  determine  if  we  can  measure 
sidelobe  levels  in  the  -60  dB  range,  and  will  modify  it  accordingly  as  needed. 
CONCLUSIONS 

1 .  Refraction  is  a  primaiy  cause  of  wavefront  distortion. 
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2. 


Refraction  seriously  distorts  the  modulus  (amplitude)  of  the  wavefront  as  well  as 
the  phase. 

3 .  Aberration  correction  algorithms  must  address  amplitude  distortion  as  well  as  phase 
distortion. 

4 .  Prior  aberration  distortion  research  has  concentrated  only  on  phase  aberration 
correction. 

5 .  Newer  and/or  stronger  algorithms  are  required  for  full  wavefront  aberration 
correction  to  achieve  the  desired  resolution  (0.2  mm)  at  the  chest  wall  and  -60  dB 
sidelobe  level. 

6 .  2-D  arrays  and  2-D  algorithms  are  required  because  wavefront  distortion  is  a  2-D 
phenomenon. 

7 .  Large  (up  to  10  cm),  highly  thinned  2-D  arrays  are  necessary  to  provide  the  desired 
lateral  resolution  and  for  practicality  and  cost. 

8 .  The  current  work  favors  the  highly  thinned,  deterministic,  multiple  ring  array  over 
the  equally  highly  thinned  random  array. 
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APPENDIX  A:  PROJECT  HIGH  RESOLUTION  ULTRASOUND  MAMMOGRAPHY 
DETAILS  OF  THE  ANALYSIS  OF  THE  HUP  DATA 

Refraction  vsScattering  Study 

Studies  in  our  laboratory  [8-9]  and  in  others[10-ll]  have  shown  that  refraction  is 
an  important  wavefront  distortion  source  inside  the  female  breast.  Refraction  is  caused  by 
speed  mismatch  across  a  tissue  interface,  e.g.,  glandular/fat  tissue.  Refraction  can  cause 
ray  bending  and  splitting.  Refracted  rays  can  be  further  refracted  at  subsequent  interfaces 
and  arrive  at  the  receiving  transducer  from  different  directions.  Such  subteams  interfere 
and  can  cause  cancellation  in  high  energy  portions  of  the  aperture  [8].  Another  wavefront 
distortion  source  is  scattering  which  is  caused  by  local  speed  variations. 

In  the  scattering  case,  the  resulting  field  E  at  the  receiver  after  a  plane  or  spherical 

wave  propagates  through  the  breast  is  the  superposition  of  the  incident  field  E-  plus  the 

scattered  field  If  the  scattered  field  is  weak,  i.e.,  E,  »  E,,  E  is  dominated  by  E-.  The 
wavefront  amplitude  distribution  is  Rician.  The  phasefront  of  E  is  primarily  the  phasefront 

of  E,  plus  perturbations  caused  by  the  scattered  field.  If  the  scattered  field  is  strong,  E- 

«E,  the  incident  field  is  lost  and  E  is  approximately  E,.  For  a  medium  containing 
multiple  scattering  bodies  or  diffuse  scatterers,  the  scattered  field  is  the  coherent  sum  of  the 
scattered  fields  radiated  from  these  scatterers.  The  central  limit  theorem  ensures  that  the 
probability  distribution  of  wavefront  amplitude  is  Rayleigh  and  phase  is  Uniform  [0,  2  tr] 
[12]. 

In  the  refraction  and  strong  scattering  cases,  the  incident  wave  may  experience 
multiple  refraction  and  encounter  multiple  scattering  bodies.  The  resulting  field  is  also 

dominated  by  E,  and  the  wavefiont  amplitude  distribution  is  also  Rayleigh. 

The  latter  case  is  what  happens  in  the  female  breast.  Wavefront  measurements  of  a 
large  population  show  that  the  wavefront  amplitude  distributions  are  approximately 
Rayleigh  at  3  and  4  MHz  (Figure  1).  Although  the  amplitude  distribution  alone  can  not 
separate  strong  scattering  from  refraction,  the  deep  nulls  frequently  observed  in  the 
wavefiont  amplitude  profiles  and  which  contribute  to  the  amplitude  distribution  close  to  the 
origin  indicate  that  refraction  is  the  likely  cause  of  this  destructive  interference. 

The  phase  aberration  compensation  algorithms  that  have  teen  developed  for  weak 
scattering  are  not  powerful  enough  when  significant  wavefront  amplitude  distortion  is 
present.  These  measurements,  therefore,  imply  that  stronger  wavefront  deaterration 
algorithms  should  be  designed  for  use  with  Rayleigh  wavefront-amplitude  distortion. 

Wavefront  Amplitude  Distortion  and  Image  Sidelobe  Levels 

Refraction  and  strong  scattering  cause  wavefiont  amplitude  distortion,  in  addition  to 
phasefront  distortion.  The  most  damaging  effect  of  severe  amplitude  distortion  is  the  rise 
of  the  sidelobe  level  in  the  image,  even  after  phase  aberration  correction  is  successfully 
applied.  High  sidelotes  can  destroy  the  desired  contrast  in  the  image  and  make  it 
impossible  to  distinguish  between  normal  tissue  and  abnormal  tissue,  e.g.,  small  tumors 
and  cysts. 

Our  in  vivo  experiment  has  shown  that  the  rise  in  the  average  sidelobe  floor  ( ASF ) 
in  a  single-source  image,  after  complete  phasefront  compensation,  is  proportional  to  the 

wavefront  amplitude  variance  cr^  normalized  to  the  square  of  its  mean  value  a  (normalized 
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variance  )  at  the  receiving  aperture  and  inversely  proportional  to  the  effective  number 
of  array  elements  N  (figure  2). 


Ampliuide 


Figure  1.  Histograms  of  wavefront  amplitude.  Dashed  lines  represent  Rayleigh 
distribution,  (a)  3  MHz,  sample  size  19013.  15  women  (30  breasts).  4  MHz, 
sample  size  2502. 4  women  (8  breasts).  (From  [3]). 

r  r 

The  equation  is  ASF  =  dB~  lOlog  ASF  =  101og(-^)  +  101og(^=^) , 

where  g  is  a  geometric  factor  relating  to  the  source  and  equals  unity  for  the  ideal  point 
source,  g  is  larger  than  unity  for  any  other  source  distribution.  The  data  were  obtained 
from  1700  samples  of  44  breasts.  The  result  is  consistent  with  the  prediction  of  the  theory 
[4].  The  measured  average  ASF  of  Fig.  2  is  -22  dB.  The  peak  sidelobe  level  is  about  7 
dB  higher  [5].  These  levels  would  be  higher  if  residual  phase  errors  remained  after 
phasefront  compensation.  The  total  number  of  array  elements  used  is  64. 
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SIDELOBE  LEVEL  vs.  WAVEFRONT  DISTORTION 


y  -  0.000529  >  0.93x 


Nomialized  Amplitude  Variance  of  Wavefront(  ^  ) 

Figure  2.  Average  sidelobe  floor  vs.  variance  of  modulus 
samples  obtained  from  44  breasts.  Transmitting  frequency  is  3MHz. 

Sidelobes  must  be  maintained  to  the  order  of  -60  dB  or 


of  wavefront.  1700 
(From  [5]) 

lower  to  distinguish 


between  small  cysts  and  tumors  of  the  same  size.  The  measured  mean  value  of  was 

(or 

0.27  or  -6  dB  and,  therefore,  it  was  shown  from  the  above  equation  that,  for  an  array  of 
100  elements,  the  expected  sidelobe  level  is  approximately  -26  dB.  There  are  two 
approaches  to  reduce  this  level  to  the  -60  dB  neighborhood.  The  first,  discussed  in  the 
next  paragraph,  is  to  increase  the  total  number  of  array  elements;  the  second  is  to  reduce 
wavefront  amplitude  distortion. 


The  relationship  between  ASF  and  was  developed  for  one-way  propagation. 

When  applied  to  echo  scanner  design,  this  model  implies  broad-beam  insonifrcadon.  This 
design  permits  the  simplest  form  of  adaptive  weight  control  because  only  the  receiving 
weights  need  to  be  controlled.  To  reduce  the  sidelobe  floor  to  a  useful  level,  the  tot^ 
number  of  array  elements  deployed  must  be  lO^.  This  would  drop  the  floor  another  30  dB. 
A  2-D  design  is  a  way  to  preserve  the  scale  of  the  transducer  while  increasing  the  total 
number  of  array  elements. 

The  second  approach  is  to  reduce  wavefront  amplitude  distortion  by  algorithmic 
design  and  transducer  shaping. 

Algorithmic  design:  It  is  possible  that  transducer  weights  of  transmission  and  reception  can 
be  iteratively  adjusted  so  that  the  transmitting  beam  pattern  and  the  receiving  beam  pattern 
are  identical.  The  result  would  be  a  squaring  of  the  side  radiation  pattern  or  a  doubling  of 
the  level  in  decibels.  Thus,  to  achieve  -60  dB  sidelobe  level,  the  sidelobe  level  of  the  array 
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need  be  no  lower  than  -30  dB.  The  iterative  phase-deaberration  procedure  developed  in 
[  1 3]  belongs  to  this  category  and  was  reported  to  be  successful  in  liver.  This  is  a  highly 
desirable  condition.  To  work  properly,  the  initial  transmitting  beam  must  have  a 
reasonably  good  pattern;  it  must  have  a  well  defined  mainlobe,  although  its  sidelobe 
radiation  pattern  can  be  high.  This  condition  is  achieved  when  wavefront  amplitude 
distortion  is  less  than  about  0.35  rms  [14].  Therefore,  in  the  case  of  a  complicated 
wavefront  resulting  from  medium  distortion,  a  satisfactory  initial  transmitting  beam  may  or 
may  not  be  formed  and,  as  a  result,  the  iterative  procedure  may  or  may  not  converge.  If 
this  procedure  fails  stronger  algorithms  which  will  correct  both  amplitude  and  phase 
distortion  need  to  be  invented. 

Transducer  shaping:  Wavefront  amplitude  distortion  can  be  reduced  through  transducer 
shaping.  Studies  in  other  laboratories  [15][10]  and  in  ours  [16]  have  shown  that  refraction 
at  skin  surface  and  at  interface  between  subcutaneous  fat  and  glandular  tissue  can  be 
reduced  significantly  by  pushing  a  transducer  against  the  breast.  But  when  a  large  planar 
transducer,  comparable  to  the  size  of  the  breast  is  employed,  it  is  difficult  to  keep  good 
contact  between  the  transducer  and  the  breast  without  extreme  compression.  The  result  is 
that  only  a  portion  of  the  transducer  is  effectively  used.  Furthermore,  extreme  compression 
distorts  the  geometric  relationships  between  tissue  types  in  the  interior  and  makes  locating 
suspicious  findings  difficult  [17].  This  is  especially  true  in  ultrasound-guided  biopsy 
where  exact  location  of  the  lesion  is  essential.  Thus  the  dual  need  is  for  a  large  2-D 
transducer  array  suitable  to  achieve  high  lateral  resolution  and  low  sidelobe  level  and  with  a 
geometry  that  minimizes  the  refractive  effect.  A  strong  candidate  is  the  cup  array  proposed 
in  [16].  The  array  consists  of  a  set  of  concentric  rings  with  rings  pulled  out  along  the 
propagation  direction  Z.  The  elements  on  each  ring  are  equally  spaced  in  angle.  This  array 
inherits  the  advantages  of  the  2-D  planar  array  and  provides  the  useful  concave  curvature  to 
minimize  refraction  effect. 

The  approach  that  aims  at  reducing  wavefront  amplitude  distortion  is  more  efficient 


than  the  one  that  increases  the  total  number  of  array  elements.  For  example,  if  -=^  were 

(a) 

reduced  by  one  order  of  magnitude,  the  total  number  of  array  elements  could  be  dropped  by 
one  order  of  magnitude  and  the  ASF  would  remain  the  same.  This  reduction  in  total 
number  of  array  elements  and  corresponding  channels  is  significant  because  the  total 
system  cost  and  technical  difficulties  of  making  such  a  large  array  will  be  reduced 
dramatically.  The  state-of-art  1-D  commercial  arrays  are  being  made  with  a  few  hundred 
elements  and  the  same  number  of  parallel  channels.  There  is  nothing  in  principle  that 
would  prevent  transducer  technology  to  extend  arrays  to  the  order  of  10^  elements  [18]. 
But  system  cost  associated  with  the  number  of  channels,  not  transducer  technology,  is  the 
limiting  factor. 

One  of  our  current  approaches  is  aimed  at  reducing  wavefront  amplitude  distonion 
through  algorithmic  design  and  transducer  shaping.  In  the  event  that  several  hundred 
elements  array  would  not  give  the  desired  sidelobe  level  with  both  efforts,  the  total  number 
of  array  elements  will  be  increased  and  optimally  deployed. 

Lateral  Resolution  as  a  Function  of  Array  Size 

In  the  absence  of  medium-induced  distortion  the  beamwidth  of  a  transducer  (in 
radians)  is  the  ratio  of  wavelength  A  to  transducer  size  L,  and  the  lateral  resolution  AS  = 
R  A  A-  is,  for  fixed  wavelength  and  depth  of  penetration  R,  set  by  the  transducer  size. 
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The  formula  AS  =  RX/L  predicts  how  lateral  resolution  should  vary  with 

transducer  size.  But  studies  in  our  laboratory  [7]  and  in  others  [  19-22]  have  shown  that  AS 
does  not  follow  the  prediction.  Figure  3  shows  how  it  actually  varies  in  in  vivo  breast. 
The  lobe  width  of  the  image  of  a  single  transducer  element  was  measured  in  a  single-source 
transmission  experiment  as  a  function  of  receiving  transducer  array  size  L.  The 
transmitting  frequency  was  3  MHz.  The  breast  was  immersed  in  a  reservoir  of  oil  whose 
acoustic  speed  match^  that  of  the  subcutaneous  fat.  Propagation  was  parallel  to  the  chest 

wall  and  2-3  cm  from  it.  The  ordinate  is  the  reciprocal  beamwidth  Au~\  labeled 
resolution.  The  curves  marked  "O"  show  the  measured  -3  and  -5  dB  lateral  resolutions  vs. 
transducer  size  for  the  oil  path;  and  similarly  for  the  curves  marked  "T"  for  resolution  in 
tissue. 

Two  breakpoints,  seen  in  Figure  3,  are  observed  in  most  tissue  curves.  Aw~‘  is 
approximately  linear  with  L  below  the  first  breakpoint;  in  that  region  the  slope  is  nearly 
constant  and  generally  smaller  than  for  the  oil  curve.  The  measured  differences  in  slope 
were  24-32%.  Beyond  the  first  breakpoint  the  resolution  flattens,  the  image  deteriorates, 
the  sidelobes  rise  and  the  lobe  representing  the  source  image  often  reduces  in  amplinide  and 
shifts  in  position.  The  first  breakpoint  is  identified  by  the  intersection  of  a  linear  fit  to  the 
left-hand  portion  of  the  data  and  a  zero  slope  fit  to  the  right-hand  portion.  Beyond  some 
aperture  size  it  is  no  longer  possible  to  unambiguously  identify  the  source  lobe.  We 
identify  this  size  as  the  size  at  which  the  peak  sidelobe  exceeds  -5  dB.  We  call  this  size  the 
second  breakpoint.  Sometimes  the  two  breakpioints  merge.  The  average  result  for  more 
than  300  cases  is  that  contrast  resolution  vanishes  beyond  L  ~  1.5  cm  and  no  identifiable 
image  exists.  This  size  is  typical  of  current  echo  scanner  transducers. 


Figure  3.  Resolution  vs.  transducer  size  L.  Transmitting  frequency  is  3  MHz.  O. 
3dB  and  0.5dB  resolution  measurements  obtained  through  oil  increase  linearly  with  L,  as 
expected.  T-3dB  and  T-5dB  are  obtained  through  breast  ID#46L.  Two  breakpoints, 
defined  in  the  text,  are  generally  observed  in  each  T  curve.  (From  [7]) 
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The  statistics  of  useful  1-D  aperture  size  provides  information  about  the  limit  for 
which  increasing  apenure  size  will  not  improve  lateral  resolution.  Above  this  limit 
wavefront  distortion  must  be  compensated. 

Wav^rorU  Correlation  Distance 

An  important  measure  is  the  correlation  coefficient  of  two  wavefronts  as  a  function 
of  lateral  distance  between  the  two  sources  [6].  It  is  shown  in  [16]  that  the  correlation 
coefficient  of  two  wavefronts  is  equal  to  the  correlation  coefficient  of  complex  images 
formed  from  the  corresponding  wavefronts.  Correlation  distance  p  is  defined  as  a 
distance  between  two  sources  when  the  wavefront  correlation  coefficient  drops  to  50%.  p 
is  a  valuable  parameter  for  it  is  related  to  how  far  one  can  compensate  wavefront  distortion 
by  weight  vector  control  at  the  aperture  and  still  form  a  good  image. 

Figure  4  shows  the  composite  autocorrelation  function  obtained  from  14 
premenopausal  fatty  breasts.  Abscissa  is  the  lateral  distance  between  the  sources  and 
ordinate  is  correlation  coefficient.  The  standard  deviations  are  shown  by  error  bars,  p  is 
about  2  mm,  which  is  exceedingly  small.  The  implication  of  the  small  correlation  distance 
is  that  weight  vector  at  the  aperture  has  to  be  changed  several  times  in  scanning  a  1cm  target 
located  at  12  cm  depth. 

One  concern  about  this  small  correlation-distance  measurement  is  the  possible 
subject  motion  which  may  have  partially  contributed  to  the  decorrelation  of  the  wavefronts. 
This  matter  will  be  evaluated  in  ongoing  experiments  using  in  vitro  measurements  (see 
Discussion). 


Figure  4.  Composite  acf  obtained  from  14  premenopausal  fatty  breasts.  Abscissa  is 
the  lateral  distance  between  two  sources  and  ordinate  is  correlation  coefficient 

Discussion 

The  fact  that  wavefront  amplitude  is  severely  distorted  is  obtained  from 
measurements  of  a  large  1-D  array  having  a  elevation  dimension  of  10  mm.  The 
propagation  path  length  was  12  cm.  In  future  3-D  echo  imaging,  the  situation  may  differ. 
First,  2-D  arrays  generally  sample  more  finely  in  elevation  so  as  to  be  able  to  form  narrow 
elevation  beams  and  thereby  resolve  elevation  multipath  and  reduce  wavefront  interference. 
Second,  propagation  path  length  (2- way  round  trip)  is  generally  less  than  12  cm  with  the 
exception  of  the  chest  region  and,  therefore,  less  severe  wavefront  distortion  might  be 
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generally  expected.  The  latter  speculation  is  possible  to  assess  by  conducting  2-D  in  \iiro 
transmission  measurements  where  the  thickness  of  the  tissue  sample  can  be  controlled. 
The  former  can  not  be  evaluated  directly  without  construction  of  a  2-D  scanner.  But  the 
characteristics  of  wavefront  interference  over  elevation  can  be  deduced  from  2-D 
transmission  measurements.  For  example,  if  scattering  were  the  dominant  distortion 
mechanism  the  wavefront  measured  at  each  elevation  would  be  statistically  the  same. 

Recently,  2-D  large-apenure  transmission  measurements  have  been  made  at  Prof. 
Robea  Waag's  lab  at  Rochester  University.  The  basic  2-D  measurement  system  was 
described  in  detail  in  [23-24].  Briefly,  wideband  3.7  MHz  ultrasonic  waves  radiated  from 
a  hemispherical  transmitter  and  the  wavefronts  after  propagating  through  breast  specimen 
were  measured  by  a  large  1-D,  horizontal  receiving  aperture.  The  array  was  translated  in 
the  elevation  dimension  to  form  a  large  2-D  aperture  (96mm  by  46mm).  The  element 
spacing  was  0.7  mm  in  the  array  dimension  and  1.4  mm  in  elevation.  Breast  samples  with 
various  thickness  were  placed  between  the  source  and  the  receiving  apenure.  Figure  5 
shows  the  measured  wavefronts  at  one  elevation  in  water  (part  a)  and  in  breast  (part  b). 
The  abscissa  is  sample  index  in  time  (I  sample  interval  corresponds  to  0.05  )is)  and  the 
ordinate  is  element  number  in  array  direction.  Near-field  curvatures  of  wavefronts  due  to 
geometry  were  removed.  Wavefront  amplimde  distortion  is  easily  observed  from  the 
discontinuities  of  the  wavefront  shown  in  Fig.Sb. 


Fig.5.  Wavefront  measurements  through  (a)  water  (b)  4  cm  thick  breast  specimen. 


Multiple  measurements  were  also  made  on  two  breast  specimens  where  the 
thickness  of  each  specimen  was  gradually  reduced  by  approximately  1  cm.  The 
progressive  distonion  in  the  measured  wavefront  with  breast  thickness  is  shown  in  fig.6. 
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One  quantitative  measure  is  the  RF  pulse  similarity  (RPS)  as  introduced  in  [25]. 
RPS  is  defined  as  the  correlation  coefficient  between  the  ideal  pulse  and  the  pulse  received 
through  tissue.  An  ideal  pulse  could  be  obtained  from  water  measurement.  High  RPS 
implies  that  the  pulse  received  through  tissue  preserves  its  shape  although  its  magnitude 
and  arrival  time  may  differ  from  what  it  should  be  in  water.  If  the  majority  of  the  RF 
pulses  received  through  tissue  has  high  RPS  and  the  distortion  is  primarily  in  arrival  time, 
time  shift  compensation  or  phasefront  compensation  may  be  enough  to  correct  for 
wavefront  distortion.  Figure  7  shows  the  percentage  of  RF  pulses  with  RPS  greater  than 
0.9  (+)  and  0.8  (*)  as  a  function  of  propagation  depth.  The  pulses  were  received  through  a 
single  breast  specimen  with  thickness  progressively  reduced  in  steps  of  approximately  1 
cm.  Extrapolations  of  the  two  curves  indicate  that  about  50%  of  the  wavefronts  would 
have  greater  than  0.8  correlation  coefficient,  and  15%  greater  than  0.9,  for  propagation 
distance  up  to  5  cm.  This  distance  is  the  order  of  the  average  penetration  distance  of  the 
breast  under  compression.  The  relationship  between  the  percentage  of  pulses  with  greater 
than  0.9  or  0.8  correlation  coefficient  and  main  beam  reconstruction  and  sidelobe  levels  of 
the  source  image  will  be  studied.  Other  quantitative  measures,  such  as  rms  arrival  time 
perturbation  vs.  thickness,  will  also  be  studied. 

This  study  will  provide  information  about  useful  penetraaon  depth  that  phasefront 
deaberration  algorithms  can  be  successfully  used  to  form  a  good  image. 
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Figure  7.  Percentage  of  RF  pulses  with  correlation  coefficient  greater  than 
0.9  (+)  and  0.8  (*). 

Summary  of  Analysis 

1.  Wavefronts  propagated  through  in  vivo  breast  and  measured  by  a  large  1-D 
aperture  are  severely  distorted.  Wavefront  amplitude  distribution  indicates  that  strong 
scattering  and  refraction  inside  breast  are  likely  causes  of  such  distonion. 

2.  Adaptive  phasefront  deaberration  algorithms  developed  so  far  can  correct 
phasefront  distortion  but  not  wavefront  amplitude  distortion,  which  causes  high  image 
sidelobes. 

3.  The  useful  aperture  size  is  limited  to  a  few  ems  or  less  in  the  breast  imaging. 

4.  The  scanning  angle  of  a  deep  sited  target  is  limited  to  few  degrees  under  weight 
vector  control  at  the  aperture. 
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